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ABSTRACT 

— 1  The  numerical  solution  of  partial  differential  equations  in  unbounded 
domains  requires  a  finite  computational  domain.  Often  one  obtains  a  finite 
domain  by  introducing  an  artificial  boundary  and  imposing  boundary  conditions 
there.  This  paper  derives  exact  boundary  conditions  at  an  artificial  boundary 
for  partial  differential  equations  in  cylinders.  An  abstract  theory  is 
developed  to  analyze  the  general  linear  problem.  Solvability  requirements  and 
estimates  of  the  solution  of  the  resulting  finite  problem  are  obtained  by  use 
of  the  notions  of  exponential  and  ordinary  dichotomies.  Useful 
representations  of  the  boundary  conditions  are  derived  using  separation  of 
variables  for  problems  with  constant  tails.  The  constant  tail  results  are 
extended  to  problems  trtiose  coefficients  obtain  limits  at  infinity  by  use  of  an 
abstract  perturbation  theory.  The  perturbation  theory  approach  is  also 
applied  to  a  class  of  nonlinear  problems .  General  asymptotic  formulas  for  the 
boundary  conditions  are  derived  and  displayed  in  detail. 
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SIGNIFICANCE  AND  EXPLANATION 


In  various  applications  tha  naad  to  solve  boundary  value  problems  for 
partial  differential  equations  in  unbounded  domains  arises.  Important 
examples  include  the  problems  of  fluid  flow  in  long  channels  and  past 
bodies.  In  order  to  solve  such  problems  numerically,  it  is  necessary  to 
compute  on  a  finite  mesh.  One  approach  to  this  problem  is  to  introduce  an 
artificial  boundary  and  solve  the  problem  on  the  finite  domain  which  it 
creates.  Here,  we  consider  the  question  of  what  boundary  conditions  to  impose 
at  the  artificial  boundary. 

First,  we  develop  an  abstract  framework  for  the  analysis  of  the  problem, 
discovering  a  useful  characterisation  of  the  exact  boundary  conditions  which 
should  be  imposed.  These  can  be  easily  represented  (and  implemented)  using 
separation  of  variables  for  problems  which  are  independent  of  the  unbounded 
coordinate.  For  other  problems,  including  a  large  class  of  nonlinear  ones,  we 
develop  a  perturbation  theory  which  leads  to  general  asymptotic  expansions  of 
the  exact  boundary  condition.  These  expansions  are  presented  in  detail. 

The  theory,  as  presented  in  this  paper,  is  mainly  applicable  to  elliptic 
boundary  value  problems  (and  ordinary  differential  equations).  However,  in 
other  works,  the  authors  extend  these  ideas  to  parabolic  and  hyperbolic 
problems .  Also,  some  numerical  experiments  will  be  presented. 


EXACT  BOUNDARY  CONDITIONS  AT  AN  ARTIFICAL  BOUNDARY 
FOR  PARTIAL  DIFFERENTIAL  EQUATIONS  IN  CYLINDERS 


Thomas  Hagstrom  and  H.  B.  Keller* 

1 .  Introduction 

Many  of  tha  boundary  valua  problems  arising  in  applied  mathematics  ara  given  on 
unboundad  doaaiM.  Examples  includa  tha  problems  of  fluid  flow  and  wava  propagation  in 
channala  or  past  bodias.  Tha  nuaer ical  solution  of  these  problems ,  however,  requires  a 
finite  domain.  In  this  paper,  wa  develop  a  theory  for  the  exact  reduction  of  a  boundary 
valua  problem  for  a  partial  differential  aquation  on  an  unboundad  cylindrical  domain  to  a 
problem  on  a  bounded  domain.  That  is,  an  "artificial"  boundary  is  introduced  and  tha 
proper  boundary  condition  to  be  imposed  there  is  derived.  In  other  works,  [8]  and  [9] ,  we 
use  our  theory  to  solve  nonlinear  problems  of  both  elliptic  and  parabolic  type. 

For  ordinary  differential  equations,  exact  reduction  theories  have  been  developed  by 
many  authorsi  deBoog  and  Weiss  [5],  Keller  and  Lentini  [11],  Jepson  and  Keller  [10]  and 
Markovich  [12].  Few  works  on  artificial  boundary  conditions  for  partial  differential 
equations,  on  the  other  hand,  have  discussed  exact  conditions.  An  exception  is  the  paper 
of  Oustafseon  and  Kreiss  [6],  where  the  form  -if  the  proper  conditions  for  a  general  hyper¬ 
bolic  problem  is  derived.  They  go  on  to  find  representations  of  tha  exact  conditions  in 
various  simple  cases  for  probleau  of  both  hyperbolic  and  elliptic  type. 

We  illustrate  the  derivation  of  exact  conditions  with  the  following  example i 

a)  V2u  ♦  a(x,^)u  -  f(x,j),  (x,^)  e  [0,»>  x  Si,  Q  c  h"  1  , 

b)  c(x'X)  <X'X>  +  «3(x,^)u(x,^)  -  Yn X  *  30  » 

c)  «<X>  <°'X>  +  b<X)u<0,x>  -  70(Z>#  X  e  n  » 

(1.1) 

d)  li*  u(x,y)  •  0  t 

X-M* 
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c(x,x>  “  4(x,;g)  *  x  *  xo  » 

«) 

f<x»X>  “  Vx,X>  “  °'  x  *  x0  ‘ 

(He  aay  that  (1.1)  haa  a  constant  tail,  due  to  condition  (l.le).)  Ha  choose  soae  point 
x  »  t  >  Xjj  as  the  location  of  the  artificial  boundary.  In  the  "tail",  x  >  T,  we  have i 

2 

a)  7  u  +  a„(X>u  "  0,  (x,X)  e  (t,»)  x  0  » 

(1.2)  b)  c,,^)  (x,^)  +  <1.(X)«<X*X>  -  0,  ^  e  3D  » 

c)  lie  u(x,g)  -  0  . 

x+“ 


Problem  (1.2)  can  be  easily  analysed  by  separation  of  variables.  Consider  the  following 
eigenvalue  problems 

.)  V*  y^,  +  .-(x,Yn(z)  -  »nYn(x),  )-ell  I 

(1.3)  b)  cm(X)  Yn(X>  +  <J„  (*)*„<*>  » 

O  /$,<%  “  1  • 

Given  certain  assumptions  on  the  boundary  condition,  (1.3b),  the  set  of  eigenfunctions, 

{Y  },  is  complete  in  that  subspace  of  L_(0)  consisting  of  functions  Mtisfying  it.  (See 
n  * 

Berezanskii  [3].)  For  simplicity,  we  further  assume  that  the  a>n  are  distinct  and  that 
»n  »  0  is  not  an  eigenvalue.  He  rewrite  the  {w^}  in  the  following  ways 

2 

w  -  o  >  0,  n  -  1,...,m  i 
n  n 

(1.4) 

2 

U)  ■  ^  Of  n  *  ^  1  f  tt+2  g  e  e  e  e 

n  n 


Expanding  u  in  tens  of  the  Yn*s, 

m 

(1.5)  U(x,jj)  -  l  cn(x)Yn(X) 

n«  1 

problem  (1.2)  becomes : 
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V  V  V 


I 

g 


v 


•fl  c  i  n  *  If « •  •  fi 

*)  c*  -  {  j  "  ,  x  e  [t,«>  > 

(1.6)  Xnon.  n  -  ^1. 

b)  11a  c  (x)  ■  0,  n  -  1,2,...  . 

x—  " 

As  (1.6a)  can  ba  trivially  aolvad,  we  aaa  that  (1.6b)  is  satisfied  it  and  only  its 

a)  c  (r)  ■  c'(t)  -  0,  n  -  1,...,a  i 

n  n 

(1.7) 

b)  c'(T)  -  -l_c_(T),  n  -  afl,  art-2,...  . 

n  n  n 

This  allows  us  to  raplaca  (1.1)  by  an  squivalant  finita  domain  problaat 

a)  V2u  ♦  a(x,2)u  -  f(x,*),  (x,*)  «  [0,t)  x  Q  , 

b)  c(x,£)  (x,jj)  ♦  d(x,x)«(*»X)  “  Y0(x,x),  X  *  8n  > 

(1.8)  c)  a(jj)  (O.jj)  +  b(x)u(0,^)  -  Y0(X>*  X  •  n  » 

d)  /ndx  u(t.X)Yb(X)  -  /jjdj  (t,x)Yn<X)  -  0,  n  -  1 . .  i 

*  ~XJa^  tt(T'*,V*)'  "  “  *f1'  “*2'*”  • 

That  is,  (1.1)  has  a  solution  if  and  only  if  (1.8)  doas  and  tha  solutions  agree  on  a  finita 
domain. 

In  section  2  of  this  work  wa  derive  boundary  conditions  for  the  reduction  of  a  general 
partial  differential  aquation  in  a  semi-infinite  cylindrical  domain  to  a  finite  one.  These 
turn  out  to  be  the  requirement  that  the  appropriate  data  at  the  artificial  boundary  lie  in 
a  certain  affine  set.  We  find  it  convenient  to  rewrite  the  problem  as  an  ordinary 
differential  equation  In  a  Banach  space,  making  transparent  the  connection  between  our 
reduction  and  the  reduction  theorems  for  the  case  of  ordinary  differential  equations.  In 
section  3  we  introduce  the  notion  of  a  dichotomy  for  our  abstract  equation  and  use  it  to 
develop  error  estimates  and  solvability  requirements  for  the  finite  problem. 

we  first  consider  the  problem  of  representing  the  boundary  conditions  in  section  4. 
Here  separation  of  variables  is  used  to  analyse  constant  tail  problems  such  as  the  one 
presented  above.  The  exact  representation  we  obtain  is  equivalent  to  (1.8d)  in  that  case. 


We  develop  a  perturbation  theory  to  analyse  non-constant  tall  probleu  In  section  5 


As  ousting  the  limiting  problea  at  Infinity  can  be  solved  by  separation  of  variables ,  a 
perturbation  expansion  of  the  exact  boundary  condition  can  be  calculated.  We  carry  out 
this  expansion  for  the  Helmholtz  equation  exterior  to  a  body,  recovering  the  conditions  of 
Beyliss,  Gunsburger  and  Turkel  [2].  Finally,  in  section  6,  nonlinear  problems  are 
considered.  We  use  the  perturbation  theory  of  the  preceding  section  to  prove,  under 
certain  conditions,  the  existence  of  an  exact  nonlinear  boundary  condition  and  to  calculate 
an  expansion  which  approximates  it. 

We  note  that  many  authors  have  derived  boundary  conditions  for  specific  problems,  we 
do  not,  in  general,  attempt  to  examine  the  connection  between  their  conditions  and  ours. 

For  more  discussion  of  these  connections  as  well  as  for  a  store  extensive  bibliography,  the 
reader  is  referred  to  Hagstrom  [7] . 


2.  Basic  Linear  Deduction  Theorem 

We  consider  abstract  boundary  value  proble 

du 


(2.1) 


l)  to  “  A(x,u  +  f<K)* 


b)  B0u(0)  -  Tf0 

c)  lim  B„u(x)  -  0 


of  the  formt 

0  <  x  <  •  » 


In  addition  we  may  imposes 

d)  lu(x)l  bounded  as  x  ♦  •  . 

For  some  Banach  space,  B,  we  seek  u(x)  e  B  for  x  e  [0,*).  we  suppose  that  A(x) ,  B0 
and  B—  are  linear  operators  with  domain  in  B,  to  which  we  also  constrain  the  range  of 
A(x).  Finally,  f(x)  6  B. 


Problems  of  form  (1.1)  can  be  obtained  from  general  partial  differential  equation 
problems  in  cylindrical  domains.  Specifically  we  consider 


(2.2) 


on  the  cylindrical  domain 


(x,x>  «  to,-)  *  0,  a  c  *n  . 

Boaogenaous  boundary  conditions  ara  imposed  on  30  involving  u>  and  its  normal 
darivativas i 


(2.3)  I  B  (2)  (X,X)  -  0,  X  e  90  . 

j-0  '3  3v3 

Ms  further  suppose  that,  subject  to  these  boundary  conditions, 

"n"  I*1 

exists  for  all  x.  Mow  (2.2)  can  be  rewritten  in  the  fore  of  (2.1a)  by  introducing 


(2.4) 


,n-2 
3  u 


3x 


n-2 


/ 

The  space,  B,  is  scsm  space  of  n-tuples  of  functions  on  0  which  satisfy  the  homogeneous 
boundary  conditions,  (2.3).  It  is  necessary  to  eliminate  inhomogeneous  conditions  on  30 
in  order  to  reduce  the  problem  to  the  abstract  form.  This  can  be  accosqslished  by  subtract¬ 
ing  a  function  that  satisfies  the  inhomogeneous  condition.  We  note  that  the  functions 

B_  affect  P-1  and,  ultimately,  A(x). 

“» j  n 

Returning  to  (2.1)  we  choose  seme  finite  point,  x  -  t,  and  attempt  to  reduce  the 
infinite  problem  on  (0,—)  to  a  finite  one  on  [0,t].  we  define  A(t »f ) ,  the  admissible 
set  of  Cauchy  data  at  x  -  t,  as  the  set  leading  to  solutions,  u,  in  the  tail, 
x  «  [t,-)«  More  precisely  we  have: 

Definition  2.5.  The  set  A(Tjf)  c  8,  the  admissible  set  at  x  -  t,  is  the  set  of  all 
Uq  «  8  such  that  there  exists  u(x)  e  B,  x  e  [t,«),  satisfying: 

a)  A(x)u  +  f(x),  t  <  x  <  —  » 


(2.5) 


b)  u(t)  -  uQ  t 


as  well  as  (2.1c,d)  as  appropriate. 
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It  is  now  possible  to  writ*  down  an  exact  reduction  of  (2.1)  to  a  problem  on  a  finite 
domain.  We  etate  the  reduction  ae  a  theorem  whose  proof  follow*  immediately  from  the 
definition  of  A(T  f  f ) . 

theorem  2.6.  Problem  (2.1)  ha*  a  solution  if  and  only  if  the  following  problem  has  a 
solutions 

dtlt 

a)  -r-  -  A(x)w  +  f(x),  x  e  [0,t)  » 

ax 

(2.6)  b)  Bq(0 ( 0 )  -  y0  , 

c)  4>(t)  e  A(x jf )  . 

Furthermore,  whenever  (2.1)  has  a  solution  u(x),  (2.6)  has  a  solution  which  is  identical 

to  u  on  [0 ,t] • 

Proofs  Suppose  (2.6)  has  a  solution,  then,  by  the  definition  of  A(T;f),  there  exists 
u+(x),  x  e  [t  ,«• ) ,  satisfying  (2.1a)  and  (2.1c,d)  as  appropriate  as  well  as  w(t)  -  U+(T). 
Define 

(«<x)  x  e  [0,T] 

u+(x)  x  e  (T,»)  . 

then,  u  is  a  solution  of  (2.1).  Now,  suppose  that  (2.1)  has  a  solution.  The  restriction 
of  u  to  [t ,*•)  satisfies  (2.5)  and  (2.1c,d)  and,  hence,  by  the  definition  of  A(T;f), 
u(t)  e  A( T ; f ) .  this  implies  that  the  restriction  of  u  to  [0,t]  satisfies  (2.6), 
completing  the  proof. 

The  set  A(rsf)  is  an  affine  subset  of  8.  A  convenient  representation  of  A  can  be 
found  in  terms  of  its  underlying  linear  subspace  and  some  particular  element  of  8 .  we 
consider  the  homogeneous  problem  in  the  tail  associated  with  (2.1): 

dv 

a)  -jj  -  A(x)v,  x  e  [T,»)  i 

b)  lim  Bmv(x)  “0  » 


(2.7) 


and,  if  (2. Id)  is  imposed 

c)  lv(x)l  bounded  as  x  ♦  • 


Corollary.  Suppose  that  for  all  vQ  e  A(t)  solutions  to  the  Cauchy  problem  defined  by 
(2.7a)  and  (2.8)  are  unique.  Then  (2.6)  has  a  unique  solution  if  and  only  if  (2.1)  does. 
Proof i  Assuming  uniqueness  of  solutions  to  (2.1)  imnediately  yields  uniqueness  for 
(2.6).  In  the  other  direction,  note  that  the  assumption  above  guarantees  the  uniqueness 
of  u+(x)  which,  combined  with  the  uniqueness  of  the  finite  interval  solution,  implies  the 
uniqueness  of  u. 

3.  Solvability  of  the  Finite  Problem 

In  this  section  tie  assume  that  solutions  to  the  homogeneous  u  problems: 

•)  ir-  A(x)v'  x0  4  x  4  xi  if  xi  >  v  xi  * x  4  x0  lf  xo  >  xi  » 

(3.1) 

b)  v(xQ)  -  v0  ; 

are  unique  for  all  XgjXj  6  [0,»).  He  define  a  solution  operator  Stx^XgiA)  in  the 
following  way: 

Definition  3.2.  Let  vQ  e  B.  If  there  exists  a  solution,  v(x)  to  problem  (3.1)  then 

(3.2)  S(x1,x0«A)v0  -  v(x1)  . 

Otherwise,  vQ  is  said  to  be  outside  the  domain  of  S(x1,x0iA). 

The  linearity  of  the  differential  equation  implies  the  linearity  of  S.  The  stated 
uniqueness  of  solutions  implies  the  consistency  of  the  definition.  Note  that  it  is 
certainly  necessary  to  restrict  the  domain  of  S  for  ill-posed  Cauchy  problems  such  as 
those  which  arise  in  the  study  of  elliptic  equations.  Whenever  S  exists,  however,  it 
does  have  the  familiar  semi-group  properties: 

•  * 

a)  S(Xj  ,x  ;A)S(x  ,xQ>  A)  «  Stx^XpjA)  ; 

(3.3) 

b)  S(Xg ,Xg|A)  -  I  . 

Hie  notion  of  dichotomies  is  very  useful  in  what  follows.  First  we  present 
definitions  of  exponential  and  ordinary  dichotomies.  These  are  adapted  from  Daletskiy  and 
Krein  [4],  with  some  modifications  required  by  the  possible  non-existence  of  solutions. 
Definition  3.4.  We  say  that  the  problem 


CWAAti 


(3.4)  ^  -  A(x)v,  x  e  [0,»)  t 

has  an  exponential  dichotomy  If,  for  any  x  e  10,-),  the  apaca  8  can  b«  decomposed  into 
a  direct  sum  of  subspaces  8_(x*)  and  8+(x*)  such  that i 

a)  If  v  e  8_(x*)  then,  for  soeie  N_  and  a_  >  0  j 

•  * 

i)  S(x,x  i A)v  exists  for  any  x  >  x  > 

(  3  e  5  ) 

#  -a_(x-x  ) 

ii)  IS(x,x  |A) vl  <  N  e  Ivl 


b)  If  v  e  8+<x  )  then,  for  some  N+  and  a+  >  0  » 

*  • 

i)  S(x,x  >A)v  exists  for  any  x  <  x  j 


-a+(x  -x) 


ii)  is(x,x  ; A) vl  <  N+e 


c)  There  exists  y  >  0,  independent  of  x  ,  such  that 

inf  „  lu  +u_l  >  y 
u±e8±(x  ) 

lu±l-1 


(This  infemum  is  typically  called  the  angular  distance  between  8+(x  )  and  8_(x  ).) 

An  ordinary  dichotoeiy  is  defined  as  above  except  that  «  0  is  allowed.  No 
"continuity"  of  the  spaces  as  functions  of  x*  has  so  far  been  required.  In  general,  we 
impose  a  sort  of  continuity  in  the  form  of  the  following  "no-mixing"  condition. 

Definition  3.8.  The  dichotomy  (3.5  -  3.7)  satisfies  the  no-mixing  condition  if  whenever 

a)  Q(x)  i8  the  projection  operator  into  8_(x) 

b)  S(Xj ,Xgj A)v  exists 


c)  Q(x1 )S(x1 ,xq(A)v  -  S(x1 ,xQ ;A)C(xq)v  . 


'  V.-.W  -  • 


%  \  V 


Assuming  thst  the  homogeneous  problem  has  a  dichotomy  in  the  tail  and  that  8_(x) 
coincides  with  the  admissible  space,  A(x),  it  is  possible  to  write  down  an  integral 


expression  for  a  particular  solution,  Up(x)  which  is  valid  whenever  lf(x)l  is 
integrable i 


(3.9)  U  (X)  -  I*  S( x ,p; A )0( p ) f ( p)dp  -  r  S(x,piA) (X-Q(p))f(p)dp  . 

P  *  " 

(The  validity  of  (2.9a)  follows  from  the  direct  differentiation  of  (3.9)  while  (2.9b)  la 
insured  by  the  identity  of  8_(x)  and  A(x)  combined  with  the  absolute  convergence  of  the 
integrals.)  Mote  that  it  is  always  the  case  that  B_(x)  •  A(x)  if  there  is  an  exponential 
dichotomy .  Then,  only  boundedness  of  If I  need  be  assuamd. 

formula  (3.9)  is  extremely  useful  in  the  developswnt  of  a  perturbation  theory,  for 
now,  we  simply  use  it  to  write  down  a  new  expression  for  the  boundary  condition,  (2.6c) » 


(3.10)  (I-Q(t))m(t)  -  -/“  S(T,p)(l-C(p))f(p)dp  . 

Extending  the  dichotomy  to  the  entire  interval,  we  now  can  prove  an  existence  theorem 
for  the  finite  boundary  value  problem  (2.6). 

Theorem  3.11.  Suppose  that  solutions  to  all  Cauchy  problems  (3.1)  are  unique  for 
x0,x1  e  [0,t]  and  that  (3.1a)  has  a  non-mixing  ordinary  dichotomy  on  (0,T)  with 
projector  Q(x)  into  8_(x).  Also  assume  that  8_<* )  *  A(t).  Then  (2.6)  has  a  solution 
for  arbitrary  f(x),  up(r)  and  yQ  in  the  range  of  B„  if  and  only  if  the  operator 


(3.11) 


♦w  -  I 


((I-Q(0))ul 
B0“  / 

has  an  inverse  with  domain  containing  all  pairs  of  the  formi 


(3.12) 


y  e  Range  (BQ) 


The  solution  is  unique  and  bounded  in  terms  of  the  inhomogeneous  data  if  and  only  if  this 
(restricted)  inverse  is. 

Proof;  We  use  the  ordinary  dichotomy  defined  by  Q(x)  to  solve  certain  initial  value 
problems.  Let 

(3.13)  W+(x)  -  S(X,T;A)(I-Q(T))up(T)  +  J*  S(x,p»A> (I-Q(p) )f (p)dp  . 

This  exists  for  all  x  on  [0,t]  by  the  definition  of  Q.  If  we  seek  solutions  to  (2.6) 


in  the  form 


(3.14) 

w(x)  ■  m+(x)  +  m_(x) 

then  w 

is  a  solution  if  and  only 

if  w_  solves 

a) 

dm 

^  -  A(x)w_  ♦  p(x)f(x)  , 

(3.15) 

b) 

B0w.(°)  -  -  *0-+(0)  » 

c) 

(l-P(T))m_(T)  -  0  . 

Me  write 

w_(x)  in  the  formi 

(3.16) 

m_(x)  •  S(x,0|A)u_(0)  ♦  Jjj  S(x,piA)Q(p)f (p)dp 

The  integral  term  again  exists  by  the  definition  of  Q  so  that  this  representation  is 
▼slid  for  any  solution  of  (3.15a).  By  (3.15c)  and  (3. Be)  we  havet 

0  -  (1-Q(t))m_(t)  -  S(t,0»A)(I-Q(0))«_(0)  » 

which,  by  the  uniqueness  of  solutions  to  the  Cauchy  problem,  implies 

( I— Q( 0 ) )w_(0)  -  0  . 

Hence,  we  can  find  a  solution  to  (3.15)  if  and  only  if  we  can  simultaneously  solve t 

(2-Q(0))tt_(0)  -  0  | 

V-<0)  myo  -  V+(0>  » 

which  in  component  form  yields  (3.11),  completing  the  proof. 

Estimates  of  the  solution  in  terms  of  the  inhomogeneous  data  are  now  obtained  from  the 
explicit  representation  in  terms  of  <fl+  and  t»_.  Assume  that 

a)  IS(x,p;A)C(p)l  «  K_(x,p) ,  0  <  p  <  x  <  T  ; 

b)  IS(x,p;A) (I-Q(p) )l  <  K+(x,p),  0  <  x  <  p  <  T  » 

(3’17)  -i/o\ 

c)  IB  III  <  K^lYl,  T  B  Range  (BQ)  » 

d)  IB0I  <  K„  . 

Then  we  have,  directly  estimating  (3.13)  and  (3.16)  and  using  ths  fact  that  Q(0)»_(0)  - 

w_(0). 
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V»  ■  m  *  -  1 


(3.18) 


)l  <  KJ*,0)K  lY  I  +  MX  lf(x)l  J*  K_(x,p)dp 
xe[0,T] 


+  MX  lf(x)l  (Jl  K  (x,p)dp  +  K  (x,0  )K  K-  /!  K  ( 0,p)dp) 
X«[0,X]  *  9 


♦  lup(X)l(lC+(x,X)  +  K_(x»0)K^K^lC+(0,X  ] 


Bquation  (3.18)  allows  us  to  sstiMta  the  arrors  caused  by  approxiMtions  to  Q(x) 


and  Up(t ) .  Suppose  we  solve  the  following  finite  problem  instead  of  (2.6): 


a)  ~  -  A(x)u  +  f(x),  0  <  x  <  t 

oX  A 


(3.19) 


b)  B0wa(0)  -  Y0  » 


C)  (I-C*(T))Mji(T)  -  (I-«*(T))U*(T)  t 


where  Q  (x)  and  u  (x)  differ  from  Q(x)  and  u  (x).  We  define  the  error,  e(x),  by 
P  P 


e(x)  =  ai(x)  -  w  (x) 


and  find  that  it  satisfies i 


a)  ^  -  A(x)a,  0  <  x  <  x 
ax 


b)  B0e(0)  -  0  , 


(3.20) 


c)  (I-Q(x) )e(x)  -  (I-Q(x))(u  (x)  -  u  (x)) 

P  P 


+  <Q(x)  -  fi*(x))(u*(x)  -  u>  (x)) 

P  * 


5  Mx) 


Wote  that  A(x),  by  construction,  is  in  the  range  of  I  -  Q(x).  (We  assusa,  of  course. 


that  a>a(x)  exists.)  Therefore  we  havei 


(I-C(x))Mx)  -  Mx)  . 


We  now  plug  into  (3.18)  to  obtain: 


(3.21) 


le(x)l  <  (K+(x,x)  +  *_(x,0)l^K0K+(0,x))IMx)l 


Further  specializing  to  the  case  of  an  exponential  dichotomy  this  becomes: 


(3.22) 


le(x)l  <  (N+e 


o+(x-x)  -a_x  -a+T 


+  H_e  K^KoN+*  )  I A  (x )  I 


That  is,  the  large  part  of  the  error  decays  exponentially  off  the  artificial  boundary. 


v 

>■ •  .  ■  .  -  -  '  r  *"  '  •  <  «*.  -  .  /,  -VA  J*.  •  •  - 


4.  Problems  with  Constant  Tails 

In  this  ssction  we  restrict  ourselves  to  problems  which  are  autonomous  in  x  for  x 
sufficiently  large.  That  is,  we  assume  there  exists  t  such  thats 

(4.1)  A(x)  5  Aw,  x  >  t  . 

We  also  require  that  the  constant  coefficient  problem  in  the  tail  be  separable.  That  is, 

we  require  that  a  complete  spectral  representation  be  associated  with 

Assumption  4.2.  There  exists  a  countable  set  of  pairs,  <*n»un) *  with  a  complex 

niwber,  Uj  f  B  and  0  not  an  accumulation  point  of  {X^}  and  there  exist  adjoint  pairs, 

(X  ,v  ),  with  v  e  Dual  (8),  satisfying 
n  n  n 

i)  A  u  -Xu  l 
*  n  n  n 

(4.2)  ii)  A  v  ■  X  v  i 

m  n  n  n 

iii>  (▼_#«_>  ”  *  . 

■  n  an 


Furthermore,  any  function  u(l  can  be  uniquely  written  in  the  forms 


(4.5)  u  -  ^  c^,  cR  -  (vft,u>  . 

Using  the  eigenfunction  expansions  defined  above,  it  is  easy  to  write  down  conditions 
for  the  existence  of  dichotomies  for  the  constant  problem  as  well  as  representations  of  the 
various  operators  discussed  in  the  preceding  sections.  In  particular  we  have  the  following 
theorem,  whose  proof  follows  imnediately  from  the  (formal)  solution  of  the  Cauchy  problem 
in  terms  of  the  eigenfunction  expansions.  (For  the  details  of  these  see  Hagstrom  [7].) 
Theorem  4.4.  a)  If  all  eigenvalues,  Xn>  of  are  bounded  away  from  the  imaginary 

axis,  then  the  homogeneous  problem  associated  with  A—  has  an  exponential  dichotomy  with 
spaces 


(4.4) 


8+  5  span{u1  s  Rt  ^  >  0}  I 
8_  S  span{u1  i  M  ^  <  0}  . 


The  exponents,  o^,  are  given  by i 
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a  -  g.l. b. |  Re  X  I  i 
RsX^>0  1 

(4.5) 

m  g.l .b. | Re  A,}  • 

R»Xt<0 

b)  Let  8+  be  defined  ee  above  and  let  8q  be  given  by: 

(4.6)  SQ  =  span{u^  i  M  •  0) 

Let  8^  •  80  be  any  direct  aim  decomposition  of  B0.  Then  an  ordinary  dichotomy  la 


induced 


by  the  apacea  8+  •  8Q  and  8_  •  8g* 


We  note  that  by  the  concluaione  of  part  (b) ,  there  can  be  many  ordinary  dichotosd.es 
aeeociated  with  a  problem  whose  operator  haa  eigenvalues  with  aero  real  part.  Which  of 
these  is  the  right  one  to  use  for  the  boundary  condition  depends  on  the  boundary  operator 
at  infinity,  B„.  Representations  of  the  solution  operator,  8,  are  also  easy  to  obtain. 

The  theorem  above  can  be  applied  to  the  example  of  section  1,  problem  (1.1). 

Rewriting  the  problem  in  first  order  form  according  to  transformation  (2.4),  the  operator 


0  -  V  -  a 
X 


Am  is  given  by: 


It's  eigenvalues  are  given  by  and  ±ia^ ,  defined  by  the  reduced  eigenvalue  problem 

(1.3)  through  equation  (1.4).  If  (1.3)  had  no  positive  eigenvalues,  the  problem  in  the 
tail  would  have  an  exponential  dichotomy.  In  the  case  of  an  ordinary  dichotomy,  the 
boundary  condition  (1.7a)  corresponds  to  the  choice: 

*0  5  V  80  5  '  • 

If,  instead  of  (l.ld),  some  other  condition  was  imposed  (for  example  a  radiation  condition) 
this  choice  would  change,  we  note  that  using  the  integral  representation  of  the  boundary 
condition,  (3.10),  the  condition  that  the  inhomogeneous  term  vanish  in  the  tail  can  be 
replaced  by  an  integrability  assumption.  Hie  boundary  condition,  (1.7),  is  then  replaced 


■x »J  KM 


by: 


which  implies 

-X  (e-T) 

(4.7)  c'(r)  “  ~X  c  (t)  -  j~  dM  f  (a),  n  »  art-1,  art-2,...  . 

n  n  n  t  n 

For  tha  imaginary  aiganvaluaa  wa  ha vat 


c;<x>] 

*  -  j  %  *•< 

-io  (a-T) 
n 

• 

1 

ian|  la  (  s-t  ) 

+  e  n 

'  1  -ia. 

c  (T  ) 

1 

4  1 

I 

A 

n 

h 

la 

n 

1  1 

ia  1 

n 

u 

which  implies 


(4.8) 


c'(t)  ■>-/**  coa[a  (s-T)]f  (s)ds  » 
n  t  d  n 

B_(T)  m  —  !“  ain[o_ (a-T )] f  (a)da  . 


a  't 
n 


For  a  ganaral  partial  differential  aquation  with  a  conatant  tail,  tha  aiganvalua 
problem  of  ita  operator,  J^,  can  be  reduced  to  an  eigenvalue  problem  for  a  partial 
diffarantial  operator.  In  particular,  it'a  aiganvaluaa,  X,  correapond  to  aolutiona  oft 


(4.9) 


(  l  Pj<X* 

j-0  3 


|j)Xj]T(X)  -  0 


coupled  with  tha  appropriate  boundary  conditiona.  Thia  ia  tha  aiganvalua  problem 
aaaociated  with  tha  Laplace  tranaform  in  x  of  the  aquation  in  the  tall,  wa  note  that  in 
practice  it  ia  tha  reduced  aiganvalua  problem,  (4.9),  which  wa  auggeat  be  eolved  to  obtain 
tha  boundary  conditiona.  The  reduction  to  firat  order  form  ia  made  in  an  effort  to 
aimplify  tha  theory.  The  uaa  of  (4.9)  to  derive  boundary  conditiona  was  firat  suggested  by 
Gustafsson  and  Kreias  [6]. 

The  completeness  of  the  eigenfunctions  of  i,  depends  on  the  completeness  of  the 
eigenfunctions  of  (4.9).  This  property  doesn't  hold  in  general  and  ia  difficult  to 
check.  For  a  class  of  elliptic  and  parabolic  problaaw,  Agmon  and  Hirenberg  [1,  Thm.  5.8] 
establish  the  completeness  of  the  eigenfunctions  and  generalised  eigenfunctions  of  (4.9) 
whose  eigenvalues  have  negative  real  part  in  the  clasa  of  aolutiona  which  are  absolutely 
lntegrable  along  with  their  first  n  -  1  x  derivatives.  In  thia  case,  the  solution  of 
(4.9)  is  guaranteed  to  yield  a  representation  of  the  admissible  space. 
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5.  Perturbation  Theory  and  Asymptotic  Boundary  Conditions 


In  th«  preceding  section  we  found  useful  representations  of  the  projection  operator, 
Q(T),  of  the  ateissible  apace  and  of  the  particular  solution,  Up(x)  for  equations  of  the 
form  (2.1)  with  constant  tails.  In  the  present  section  we  relax  this  assumption  and 
replace  it  with: 


(5.1) 

Bquivalantly  we  write t 


lim  A(x)  ■  Aw  . 
x** 


(5.2) 


A(x)  -  *  B(X)  ) 

lim  I B( x) I  -  0  . 

x*~ 


Assuming  A_  has  a  dichotomy,  it  is  possible  to  make  an  asymptotic  analysis  of  the 
perturbed  problem  defined  by  A(x).  In  particular,  we  obtain  representations  of  the 
projector,  Q(t),  into  the  admissible  apace.  Consider  the  homogeneous  problem  in  the 
tail: 

dv 

a)  —  -  J^v  ♦  B(x)v,  x  >  t  i 
(5.3)  b)  lim  B„v(x)  -  0  t 

c)  lu(x)l  bounded  as  x  ♦  •  . 


Treating  B(x)v  as  an  inhomogeneous  term,  we  have,  by  (3.10),  that  v(x)  must  satisfy: 
(5.4)  (I-q„(T))v(T)  -  -  /*  S(T,piA-)(I-^.(p))B(p)v(p)dp  . 

Also,  from  (3.10),  we  have  a  representation  of  v  which  must  be  valid  if  v  exists) 


(5.5) 


v(x)  -  S(x,TfAM)paa(T)v(T)  ♦  J*  S(x,p:A<,)Qt>(p)B(p)v(p)dp 
-  S(x,p)A-)(I-91>(p))B(p)v(p)dp  . 


Let  any  B  /^,(t)  be  given  and  replace  Q„(t)v{t)  in  (5.5)  by  CQ.  If  the  following 
condition  holds: 

(5.6)  sup  I [J*  8(x,p> A  )Q  (p)B(p)  •  dp  -  /"  S(x,p» A^) (I-Q^tp) )B(p)  •  dp) I  -  K  <  1  ; 

x>T 

then  the  contraction  ampping  theorem  can  be  used  to  establish  the  existence  of  a  unique 


bounded  solution  to  equation  (5.5), 

qjr<T,e0>  -  c0  » 

(5.7) 

<T-0„(T))v<tiCo)  - 


v(siCg).  Furthermore,  we  clearly  have  that! 


-  J"  B(T,p»\,)(I-^,(p) )B(p)v(p»C0)dp  . 


Banco,  whenever  (5.6)  la  valid,  we  can  find,  for  any  C  A-(t),  a  unique  element, 
v(t >Cq ) ,  of  A(t).  A  projector  into  A(t)  ia  given  lnplicltly  by  (5.7)t 

(5.8)  Q(T)C  -  Q„(T)C  -  /"  ^pS(T,p|Aw)(X-Q|t(p))B(p)v(p!^a(T)C)  . 

These  conditions  lead  us  to  the  following  theorem! 

Theorem  5.9.  Me  auppoae  that  either  the  unperturbed  problem  has  an  ordinary  dichotomy  and 
lB(x)l  is  lntegrable  or  that  the  unperturbed  problem  has  an  exponential  dichotomy.  Then, 
for  T  sufficiently  large,  a  unique  solution,  v(xi£q),  exists  for  any  (Q  •A-(x)  and 

(5.8)  is  valid. 

Prooft  Xt  is  only  necessary  to  satisfy  (5.6).  Xn  the  first  case  we  havei 

K  <  (W+  ♦  B_)  IB(x)ldx  | 

while  in  the  second  we  havei 

B  B 

K  <  (—  +  -^)  max  IB(x)l  . 

-  *>t 

For  both  cases,  the  assumptions  on  B  allow  us  to  make  the  right-hand  sides  arbitrarily 
small  by  choosing  t  sufficiently  large,  completing  the  proof. 

The  contraction  mapping  solution  of  equation  (5.5)  leads  to  a  natural  iterative 
procedure  for  the  approximation  of  v(x>CQ)  and,  ultimately,  of  the  operator  Q.  we  leti 

v*0,(xiC0)  ■  8(x,TiA-)f(>  i 

(5.9)  v<n+1,(x»C0)  -  v(0,(x»e0>  +  J*  Ap  8(x,p»A-)Qt.(p)B(p)v(n)(p,£0) 

-  ^  <»P  8(x,piA<>)<I-91,(p))B(p)v<n,(p,C0)  . 


Then,  by  the  contraction  estimates! 


<5.10) 


i  \  _n+1  fn, 

lvln,<x,C0>  -  v<xiC0)l  <  lv(O,<x»C0)l  . 

M  define  our  nth  approximation  to  A(t ) ,  Q*n*(x),  by: 

(5.11)  S»<n)(T)^  -  QjTjC  -  /"  dp  S<T.p,A-HI-QJp))B<p>V(,,-1><p,git(T)S)  . 

The  error  due  to  thie  approximation  is  estimated  by: 


<5.12) 


IQ(t>€  -  C<n><T)C«  <  IV(0)<x»<^<t)C)I 


(Mote:  in  all  eases  the  norm  of  a  8 -valued  function  of  x  is  taken  to  be  the  maximum  in 
x  of  its  B  norms.) 

He  now  apply  these  results  to  the  case  when  the  constant  tail  problem  has  an 
exponential  dichotomy  and  A—  has  a  complete  spectrum.  He  assume  that  B(x)  has  an 


expansion  of  the  form: 
(5.13) 


.  1  _<1)  .  1  (2) 
B(x)  -  —  B  +  -T  B  +• 

x‘ 


(The  expansions  could  easily  be  carried  out  for  more  general  forms.)  Plugging  into  the 
formulas  above  we  have: 

M.  .  X  (x-t)  w  X  (x-p)  X  (p-t) 

V  ’(xiC)  -  l  cne  *  l  l  f*  dp  e  Bnm<p>cm* 

n  n  m 


(5.14) 


ReX  <0 
n 


ReX  <0  ReX  <0 

n  m 


V  v  r  Xn(x"p>  ^.(P-T) 

).  ).  <*P  •  Bnm(p>cm*  » 

n  m 

ReX  >0  ReX  <0 

n  ■ 


<vn,c>  , 


(5.15) 


B^Jx)  “  <*„#*<*)«*> 


Using  (5.13)  and  approximating  the  integrals  using  integration  by  parts  yields  to  within  an 
0(— _)  error: 


X_(X-T  ) 


-  -  -  «  .  X _  - 


(5.16) 


<x,5)  -  l 


X(X-T ) 

u  c  •  n  ♦ 

n  n 


.M) 


RaX  <0 
n 


RbX  <0  RbX  <0 

n  ■ 

VXn 


V  rm  /• 

i  un  r^r  v- 

*  n 


X  (X-T) 
n 


v  _(1)  "n'~  ' ,x.  .  r  r  me 

+  ^  V*  Cn*  1°9<t>  +  l  i  \  T-=r  c.  ~ 


X_(x-r ) 


.M) 


X_  ( X-T  ) 


RaX  <0 
n 


RaX  >0  RaX  <0 

n  a 


Putting  this  axpraaalon  into  (5.11)  and  approximating  tha  intagrala  in  a  similar  fashion 
yialdsi 

Q(T)C  -  J  e  u  ♦  l  X«  (B*1*  +  -  B(2,)c  1  --  - 

L  nn  L  “  n  n>  t  nm  m  (X  -X  )t 


(5.17) 


RaX_  <0 
XI 


R»X  >0  RaX  <0 
n  bi 


-  I  l  «-(1> 


.  _  c 

n  ni  ■ 


RaX  >0  RaX  <0 

n  m 


2  2 

<WT 


-  y  y  y  «  b(1,b(1,c 

t  J  t  Vn)  jm  a 


RaX  >0  RaX  <0  RaX  <0 
n  3  a 


nj  1*  T2(X  -X  )(X  -x  ) 
■  n  3  n 


♦  I  l  l  «„*il>81m>cm  - - - ♦  O(-lj)  i 

RaX  >0  RaX  >0  RaX  <0  3 

11  J  » 


en  *  (V5)  ‘ 

Tha  ganarality  of  tha  axpansion  givan  above  aakas  its  automatic  computation  a  raal 
possibility.  Mota  that  tha  axpansion  is  equivalent  to  tha  ona  obtain ad  by  Japson  and 
Keller  [10]  for  ordinary  diffarantial  aquations. 

Formula  (5.17)  can  ba  appliad  to  tha  laplacian  example ,  (1.1),  where  tha  potential 
a(y)  is  raplacad  byi 

(5.18)  a(x,j[)  “  +  ■;  *!<*>  ♦  *2(X>  ♦•••  • 

x* 

(i) 

Than,  tha  matrix  alaaants  ■  are  given  byi 


(5.19) 


B«  mjirfo*x  *a<z>*a<z>\W  • 

Ik 

BXpansions  of  a  particular  solution  can  be  derived  in  a  similar  manner,  let  t^(x) 
be  any  particular  solution  of  the  unperturbed  problem.  Then,  a  solution  of  the  integral 
aquation i 

(5.20)  u(x)  -  u„(x)  +  J“  S(x,p)Q„(p)B(p)u(p)dp  -  /“  S(x,p)  (I-Q^(p)  )B(p)u(p)dp 
is  a  particular  solution  of  the  perturbed  problem.  Given  the  inequality  (5.6),  a  unique 
bounded  solution  of  (5.20)  exists  by  the  contraction  mapping  theorem.  It  can  be 
approximated  by  an  iterative  process  analogous  to  the  one  described  by  equation  (5.9). 
Perturbations  of  the  inluxaogeneous  term  could  also  be  included. 

Finally,  we  note  that  (5.17)  is  valid  for  some  problems  which  do  not  satisfy  (5.6). 

An  Important  example  is  afforded  by  the  exterior  Helmholtz  problem  in  two  dimensions.  The 
equation  in  the  tail  iss 


(5.21) 


2  2 
3  u  1  3u  1  3  u  2 

2  r?7+J2iT+,tu*0,  r€  t1'"’'  9  ®  [0'210 

3r  r  36 


together  with  boundary  conditions 


(5.22) 


a)  u  periodic  in  9  i 

b)  lim  r/2  (|^  -  iku)  -  0 


Rewritten  in  first  order  form  these  become  •. 


,  *  .  ,  0  k  .1,1  0.,»,  .  1 

a)  (u)  +  (_t  Q  )(u)  +  —  <0  )(  )  +  — j 


36* 

2  0  0 


O’O  » 


(5.23) 


b)  (“)  periodic  in  6  ; 


c)  lim  r  2 (u  -  iku) 


There  are  two  obstacles  to  the  application  of  the  preceding  theory  to  problem 


(5.23).  The  first  is  that  the  perturbation 


JL  39 

_2  \  0  0 


)  ill  0  (ii 

h  v ♦  o  >0 


is  apparently  unbounded •  At*  second  is  that  the  perturbation 

if1  0) 

r  '■0  O' 

is  non-integrable  while  the  limiting  aquation , 

(5.24) 

has  an  ordinary  rather  than  an  exponential  dichotomy.  Nonetheless,  it  is  possible  to  apply 
formula  (5.17),  or  any  higher  order  approximation  to  the  boundary  condition,  to  this 
problem.  (It  is  necessary,  of  course,  to  identify  the  imaginary  eigenvalue  ik  with 
eigenvalues  with  negative  real  part  and  -ik  with  eigenvalues  with  positive  real  part  when 
applying  the  formulas.)  The  resulting  boundary  condition  isi 


(5.25) 
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2  2 
4kV 


(W(T) 
U<T  ) 


)  .  J_  (1  ♦  _L_ 

1  2ikr  v2  2ikt  ae2;|  1 


1  -lk\ 

2  2 

2 

2 


2ik 


(“,T>) 
u(T ) ' 


(“<T>) 
u(T  )  ' 


M  (  T  )  1  0 

1  1  |®  +0t7»  -  ‘o’ 


(W(T)) 

'u(T)' 


,  2ik 


which  can  be  written i 
3u 


(5.26) 


1 


{ t ,0 )  -  lku( t ,0 )  -  u(t,0)  - 


1  »2u  ,  1 

—  — 2  (T'0) - 2 


U(T,0) 


2ikT  30  8ikr" 

The  validity  of  (5.26)  can  be  established  by  other  means.  See,  for  example,  Bayliss, 
Guncburger  and  Turkel  [2].  We  note  that  the  error  depends  on  higher  0  derivatives  of  u. 


6.  Nonlinear  Problems 

In  this  final  section  we  apply  the  perturbation  theory  to  nonlinear  problems.  We 
restrict  ourselves  to  abstract  problems  of  the  formi 
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a)  —  -  F(u),  x  >  t  I 


b)  lin  u(x)  *  v^,  > 


c)  Flu.)  -  0  i 


where  u(x)  is  an  element  of  some  Banach  space,  8,  end  F  Is  a  nonlinear  operator  with 


dosMln  and  range  In  B.  Letting  v  -  u-u— ,  we  rewrite  (6.1): 


«)  -  Fu<u„>v  +  R(v),  x  >  t 


b)  11a  v(x)  -  0  i 


c)  R( v)  =  F(uw+v)  -  Fu(»^.)v 


One  approach  to  the  solution  of  (6.1)  or  (6.2)  would  be  Newton's  method.  Then,  the 


theory  of  the  preceding  sections  could  be  applied  at  each  stage  of  the  Iteration,  we. 


however,  choose  to  work  directly  with  (6.2),  deriving  exact  boundary  conditions  which  can 


be  approxiaated  by  the  methods  of  section  5. 


We  generalise  the  notion  of  an  a<hxisslble  set  (Definition  2.5)  to  be  applicable  to 


(6.2).  Note  that  it  is  no  longer  an  affine  subset  of  BL  Central  to  our  analysis  is  the 


behavior  of  solutions  to  the  linearised  problea  in  the  tail: 


d?  -  ru(u>,  x>t 


Treating  the  nonlinearity,  R(v),  as  an  inhomogeneous  term  leads  to  the  following 


equations  for  v,  which  are  analagous  to  (5.4)  and  (5.5); 


(I-Q.(t))v(t)  -  -  /“  S(T,p;Fu(u-))(I-(^t(p))R(v(p))dp  ; 

v(x)  -  S(x,T»Fu(u-)  )C>.(T)v(T)  +  J*  S(x,p;Fu(u-)  )Q„(p)R(v(p)  )dp 


-  J“  S(x,p;Fu(u-))(I-(^t(p))R(v(p))dp 


Here,  Q*  projects  into  the  admissible  space  of  the  linearised  problea  (6.3).  As  in  the 


linear  case,  the  condition  that  (6.4)  and  (6.5)  be  simultaneously  solvable  is  viewed  as  a 


condition  for  the  admissibility  of  v(t). 


«*. 


t* 


Following  the  derivation  for  the  linear  problem,  we  let  «  A.(t)  be  given  and  use 
a  contraction  argument  to  establish  the  existence  of  a  solution  to  the  integral  equation, 

(6.5) ,  with  Qaa(T)v(t)  replaced  by  £q>  tXie  to  the  nonlinearity,  some  additional 
assumptions  are  neededt 

Assumption  6.6. 

a)  There  exists  6  >  0  such  that  if  u^Uj  e  B  and  lu^l  <  8,  i  ■  1,2,  then 

sup  lj*  S(x,p»Fu(u<>)  JQ^p)  (R(u1 )  -  R(Uj  )  )dp  -  /^S(x,p»F  (u^) )  (I-Q^p) )  (RfUj)  -  R(u2)>dpl 
x>T  a  x 

<  Klu1  -  U-l,  K  <  1  . 

b)  There  exists  8^  >  0  such  that  if  u  8  B  and  lul  <8,  then 

sup  if*  S(x,p;F  (u^JlQ^lplRtuJdu  -  f“  S(x,p»Fu(u-))(I-(^(p))R(u)dpl  <  8-81  . 

x>T  U  iJ 

c)  sup  IS(X,TjFu(uto))50l  <  S,  • 

X>T 

Given  these,  a  solution  to  (6.5)  is  guaranteed  by  the  contraction  mapping  theorem. 
Denoting  this  solution  by  v(xitQ),  an  exact  boundary  condition,  valid  for  small  boundary 
data,  can  be  written  down  from  (6.4): 

(6.6)  (1-<Z'(t))v{t)  -  -  S(T,p>Fu<uw))(l-^i(p))R(v(p)q|t(T)v(T)))dp  . 

An  approximation  to  (6.6)  can  be  obtained  from  an  iterative  approximation  to  the 
solution  of  (6.5): 

a)  v<0>(x;£0)  -  S(x,T  jF^u^ )  )?Q 

(6.7)  b)  v(n+1)(x,C0)  -  v<0,(xj50)  +  J*  S(x,p,Fu(ut.))Qo,(p)R(v(n)(p;C0))dp 

-  /"  S(x,p»Fu(ut,))(I-^1>(p))R(v(n,(p»eo))dp  . 

The  nth  approximation  to  the  boundary  condition  is,  then,  given  by: 

(6.8)  (I-Q„(T))v(T)  -  -  /"  dpS(T,p»Fu(u-))(l-Q-(p))R(v(n)(p»Qc>(T)v(T)))  . 

Error  estimates  follow  as  in  the  linear  case  and  will  be  proportional  to  Kn+1lv*°*l 

-IJ 

which,  in  turn,  we  expect  to  be  proportional  to  lv(T)l  .  Note  that  R  will  often  be 


given  as  an  expansion: 


R(  v)  ~  ~  F  (u  )w  +  |f  (u)vw  +. 
2  uu  «  6  uuu  • 


'A .v.y.vy.y.v. 


He  taka  ••  many  tans  in  this  expansion  whan  evaluating  the  integrals  as  is  consistent  with 
the  number  of  tens  in  (6.7)  we  intend  to  retain. 

Assume  now  that  the  linearised  operator,  ru(ua)>  has  a  complete  spectrum.  Then,  in 
order  to  satisfy  part  (a)  of  assuaption  (6.6),  it  is  necessary  to  assume  that  there  is  an 
exponential  dichotoay.  From  (6.9)  we  derive  the  following  representation  of  R(v)  in 
tens  of  the  eigenfunctions  of  Fu(u„)  i 
m 

a)  v  "  l  cnun'  cn  “  <vn'u)  » 


(6.10) 


b)  R(v) 


l 

n-1 


Y.(v)u  , 
n  n 


Vv)  ~  l 

i,j 


a(n,c  c  + 

13  13 


l 

i,3,k 


„(n) 

6ijkcicjV“ 


c)  -S’  "  (vn'  2  WV'  "  (V  ?  Fuuuuiu jV‘ ’ * 


.(D, 


The  function  v  (x;£)  is  given  byi 

X  (x-T) 


.(1», 


V  <X|C)  -  l  u  c  e 
n 


V  .  v.  ,  in.  •  n  . 

+  l  l  l  UnBij  (l>X>f  7 


(X/Xj(x-t)  Xn(x-r) 
-  e 


i  j  n 


ReX  <0 
n 


ReX  <0  ReX. <0  ReX,<0 
n  i  j 


W*n 


X  (x~t) 

*  l  l  l  Viicici*  (X’T> 

n  i  j  n  X3  x  3 


(6.11) 


J 

ReX  <0  ReX.  <0  ReX  <0 
n  i  j 

X .  +x  «x 

i  j  n 


'III 

ReX  >0  ReX.  <0  ReX.<0 
n  l  3 


,  (X  +X.)(x-t)  , 

1 -  13  +  0(iei  )  , 


c  “  (v  .£> 
n  n 


This  yields  the  following  approximation  to  the  boundary  condition,  which  we  write  in  tens 
of  the  expansion  coefficients.  Here,  n  is  such  that  Re  X  >  0. 


(6.12) 


■  l  l  1  \  +v hr  *  l  l  £  eijkcicjck  x  Vx  lx 

i  J  ■’  i  J  n  i  j  k  313  Ai  j  Ak  n 

R*Xi<0  R«X.<0  ReX1<0  ReX j<0  ReXk<0 

l  |  I  ^aijaki  +  ajiawJ  clckci  (X.+X.-t-X  -X  XX.+X.-X  ) 

1  j  K  *  K  A  1  A  1  j  D 

R«X1<0  RaX j  <0  RaXk<0  ReXt<0 

+  i  J  k  *  +  aJi“w}  W*  <  VWV'WY  ' 

R®Xi<0  ReX^>0  RaXk<0  ReXt<0 

All  ganeral  formula  can  ba  applied,  for  example,  to  nonlinear  elliptic  probleme  of 
the  form 1 


a) 

v2u  -  f(u,^)»  {*,%)  e  tt,**)  *  o  t 

(6.13) 

b) 

Bjju  «  0,  y  e  30  » 

c) 

lim  u(x,y)  «  u,(y)  j 

where  u-(^)  satiefiesi 

(6.14) 


a)  Vjtt.  -  ftu.,!),  X  €  1)  , 

b)  -  0,  x  *  80  • 


The  linearized  equation  in  the  tail  ie  given  by  1 

(6.15)  V2v  -  fu(u„,X)v  -  0  1 

which  ia  of  the  form  analyzed  in  aection  3.  The  condition  that  (6.15)  have  an  exponential 
dichotomy  is  that  all  eigenvalues,  aR,  of  the  problem 

a)  V2Y  -  f  <u  ,v)Y  -  a  Y  ,  y  e  0  » 

^n  u  *  *  n  n  n  * 

(6.16) 


b)  BjjYn  -  0,  X  e  30  t 


be  negative.  Then,  the  following  boundary  condition  can  be  derived  from  (6.12): 


(6.17) 


v  v  n  1  j  ~  v  7  -n  ijk 

°n  “  ”\.Cn  ”  1  ^  1+1+1  "  ^  ^  ^  ^Hk  1  +1  +1  +A 

n  n  n  1-1  j-1  iJ  *i  *j+*n  i-1  j-1  k-1  13k  *i+*j+*k+*n 

r  7  *  v  n  j  cickci  r  1  1  ^ 

.  ,  A  A  A**1^  **  X. +X.+A.+A  U.+V.+X.  "  A.+A  +X.  J  ’ 

i-1  j-1  k-1  t-1  J  k  t  i  n  k  l  j  1  n  j 


n  ■  1*2,3 1« 


Here  we  have: 


J(l  <&V*,V<T'*)  » 


-°1  » 


(6.18) 


-  /a  ^  T  V*1  » 

Bijk  -  fa  *  i  fuuu(tt-'*,Vz>V*,V*,V*,  • 

The  quadratic  approximation  to  thia  condition  ia  uaed  in  a  numerical  computation  by  the 
authors  in  [8] . 
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